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Abstract. Electric-field induced magnetization switching in multiferroic magneto¬ 
electric devices is promising for beyond Moore’s law computing. We show here that 
interface-coupled multiferroic heterostructures, i.e., a ferroelectric layer coupled with 
a ferromagnetic layer, are particularly suitable for highly-dense, non-volatile, and 
ultra-low-energy computing. By solving stochastic Landau-Lifshitz-Gilbert equation 
of magnetization dynamics in the presence of room-temperature thermal fluctuations, 
we demonstrate that error-resilient switching of magnetization is possible in sub¬ 
nanosecond delay while expending a minuscule amount of energy of attojoule. 
Such devices can be operated by drawing energy from the environment without the 
need for an external battery. 
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1. Introduction 

Electric field induced magnetization switching mechanism holds profound promise for 
computing in beyond Moore’s law era in El- In multiferroic magnetoelectrics, application 
of an electric held can rotate the magnetization via converse magnetoelectric effect, 
however, such materials in single-phase usually have issues of weak coupling between 
polarization and magnetization and they operate only at low temperatures IS i. 
Although, new concepts may come along on switching in single-phase materials [51 E], 
strain-mediated multiferroic heterostructures la El [3 i consisting of a piezoelectric 
layer coupled to a magnetostrictive nanomagnet, is shown to be very effective. With 
appropriate choice of materials, when a voltage of few millivolts is applied across 
such heterostructures, the piezoelectric layer gets strained and the strain is elastically 
transferred to the magnetostrictive layer rotating its magnetization. Such switching 
mechanism dissipates a minuscule amount of energy of ~1 attojoule (aJ) in sub¬ 
nanosecond switching delay at room-temperature [9]. This study has opened up a new 
held called stmintronics miiniiii] and experimental ehorts to demonstrating electric- 
held induced magnetization switching are considerably emerging HailSllIllIH]. In a 
strain-mediated multiferroic heterostructure, the strain transferred by the piezoelectric 
layer to the magnetostrictive nano magnet, can only rotate the magnetization 90° 
in steady-state consideration [16]. Although there are proposals of 90° switching 
mechanism [m la Eg, it is explained that a complete 180° switching is possible if 
we consider the dynamics of magnetization into account rather than assuming steady- 
state scenario; basically the magnetization’s excursion out of magnet’s plane provides 
an equivalent asymmetry to cause switching in the correct direction mm- 

Although the aforesaid strain-mediated mechanism in multiferroic heterostructures 
is promising, it would be of substantial interest if there exists a strong coupling between 
the polarization and magnetization at the heterostructure interface. Recently, interface- 
coupled multiferroics are proposed based on density functional theory (DFT) of hrst- 
principles calculations [18]. Despite lacking experimental verihcation for this specific 
case as of now, considering that hrst-principles calculations have been proved to be 
very useful [19] and with the experimental progress on similar front [201 EB 1221 l23l [23] 
(also using ferromagnetic oxides [25l rather than ferromagnetic metals), there is a 
considerable interest on such coupling mechanism [THl 123 EHl EH] EHl El]- Figure El 
depicts the interface coupling between polarization and magnetization in a multiferroic 
heterostructure. The polarization direction in the P-layer determines the ground state of 
the trilayer Mi/spacer/M 2 . For polarization P|, parallel alignment (P-alignment) in the 
trilayer is preferred, while for the polarization P-i-, antiparallel alignment (AP-alignment) 
in the trilayer is preferred. This unique coupling phenomenon alongwith electric-held 
induced polarization switching makes the switching of magnetization in the Mi-layer 
non-volatile. Also, if a voltage with certain polarity is applied and maintained, the 
state of the system remains unaltered. This is advantageous over the strain-mediated 
switching, which just toggles the magnetization states. There are other exchanged 
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Figure 1. Schematics of the interface-coupled multiferroic magnetoelectric devices 
(see Ref. [IB])- The unique coupling between the polarization in the P-layer and 
the trilayer Mi/spacer/M 2 allows the polarization direction to dictate the magnetic 
ground state in the trilayer. If the polarization points downward (P^), P-alignment in 
the trilayer is preferred while an upward polarization (P.|-) prefers the AP-alignment. 
Application of a voltage with correct polarity can switch the polarization and hence 
the magnetization Mi gets switched too due to interface coupling. At the bottom 
of the figure, the axis assignment for the dynamical motion of magnetization Mi in 
standard spherical coordinate system is shown. 


coupled systems with an insulating spacer layer, but the interlayer exchange coupling 
energy is small [3^133] . Also, there are other schemes with non-magnetic spacer layer to 
preserve large interlayer exchange coupling, however, the electric held required is high 
and also the switching is volatile |34] . 

Here, we study the magnetization dynamics in the interface-coupled multiferroic 
heterostructures by solving stochastic Landau-Lifshitz-Gilbert equation in the presence 
of room-temperature thermal huctuations. Such phenomenological study of switching 
has been very useful to understand the performance metrics of magnetic devices |2l 
U I]- First, we model the interfacial anisotropy in the interface-coupled multiferroic 
heterostructures and then we analyze the dynamics of magnetization, which shows that 
switching in sub-nanosecond delay is possible while expending only ~1 aJ of energy at 
room-temperature. The strong interface anisotropy makes the switching error-resilient 
and fast, and it allows us to work with nanomagnets with small dimensions, i.e., the 
magnetization is stable with ~10 nm lateral dimensions even in the presence of room- 
temperature thermal fluctuations. Such superior performance metrics of area, delay, 
and energy are particularly suitable for computing in beyond Moore’s law era. 
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2. Model 

Figure [H shows the schematics of the interface-coupled multiferroic heterostructure 
devices and the axis assignment for the orientation of magnetization Mi. The standard 
spherical coordinate system with 6 as polar angle and 0 as azimuthal angle is utilized. 
The magnetization Mi orients along 6 = 180° if polarization points downward {Pi), 
while Ml is along 0 = 0° if polarization points upward {Pi). The elliptical cross-section 
of Ml lies on the y-z plane (0 = ±90°) with its major axis pointing to 2 ;-direction 
and minor axis in ^/-direction. The dimensions of the nanomagnet Mi along the z-,y-, 
and a:-axis are a, b, and /, respectively. So the nanomagnet’s volume = {n/4)abl. 
Magnetization’s ground states reside on the magnet’s plane {y-z plane, 0 = ±90°), 
however, during the dynamical motion, magnetization can deflect out of magnet’s plane 
and any deflection from 0 = ±90° is termed as out-of-plane excursion. 

The interface anisotropy energy in the nanomagnet Mi per unit volume is modeled 
as 

Ei{e,t) =-MHj{t)cose, (1) 

where M = fioMg, Mg is the saturation magnetization, and Hj is the interfacial 
anisotropy held. If Hj = —Hj^rnax, the ground state of magnetization Mi points 
along 9 = 180° and if we vary Hj from —Hj rnax to Hj „iax, the ground state orients 
along 6 = 0°. The total anisotropy of the magnet is the sum of the interface 
anisotropy alongwith the other anisotropies like magnetocrystalline anisotropy and 
shape anisotropy ISli, however, due to strong interfacial anisotropy compared to the 
other anisotropies, we consider only the interfacial anisotropy (i.e., Etotai — Ej) for 
brevity. We assume that the nanomagnet has a shape of an elliptical cylinder with the 
ellipse’s major axis along the z-direction, so that the magnetic easy axis becomes along 
the ±z-direction. 

The magnetization M of the single-domain nanomagnet Mi (having constant 
magnitude of magnetization but a variable direction) can be represented by the unit 
vector in the radial direction hr in spherical coordinate system {r,6,(j)), i.e., rim = 
M/|M| = hr. The other two unit vectors in the spherical coordinate system are 
and 60 for 6- and 0- rotations, respectively. The torque Tj acting on the magnetization 
due to interface anisotropy can be derived from the gradient of the energy and is given 
by 

Ti(6',t) = —rim X VEi{9,t) = —MHj{t)sin9e^. (2) 

Note that the torque Ti acts along the out-of-plane direction, so that the magnetization 
can delflect out of magnet’s plane (i.e., 0 can deflect from ±90°). 

The effect of random thermal fluctuations is incorporated via a random magnetic 
held h(t) = hx{t)ey_-\-hy{t)ey-\-hz{t)ez, where hi{t) {i = x, y, z) are the three components 
of the random thermal held in Cartesian coordinates. We assume the properties of the 
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random field h(t) as described in Ref. [3S]. The random thermal held can be written 
as [35] 


hi{t) 


2akT 

\-f\MnAt 


G(o,i)(t) (iex,y,z) 


(3) 


where a is the dimensionless phenomenological Gilbert damping parameter, 7 is the 
gyromagnetic ratio for electrons, 1/At is the attempt frequency of thermal noise, G is 
the volume, and the quantity G(o,i)(t) is a Gaussian distribution with zero mean and 
unit variance. 

The thermal held and the corresponding torque acting on the magnetization 
per unit volume can be written as Hth( 6',0, t) = Pe{0, (j),t) Bq + P^{9, (j),t) e^i, and 
Tth( 6 ', 0 , t) = rim X Hth( 6 ', 0 , t), respectively, where 


Pe( 6 ', 0, t) = M[hx{t) cos9 coscj) + hy{t) cos9sin(j) — hz{t) sin9], (4) 

P^{9,(j),t) = M[hy{t) coscj) — hxit) sincj]. (5) 


The magnetization dynamics under the action of the torques Tj and Tth 
described by the stochastic Landau-Lifshitz-Gilbert (LLG) equation as follows. 


IS 


C^llir 

dt 


a n 


X 


dt 


—-[Ti TxhI 

M ^ rnj 


( 6 ) 


Solving the above equation analytically, we get the following coupled equations of 
magnetization dynamics for 9 and (j)\ 

(1 + ^ [-aMHi{t) sin9 + {aPe{9, 0, t) + P^{9, cj), t))], (7) 

(1 + ^ M [MHj{t) - [sin9]-^ {Pe{9, 0, t) - aP^{9, 0, t))] 

{sin9 7 ^ 0). ( 8 ) 


We solve the above two coupled equations numerically to track the trajectory of 
magnetization over time, in the presence of room-temperature thermal huctuations. 

From Eqs. ([7]) and (0), we see that the torque acting in the 0-direction is much more 
higher than the torque exerted in the ^-direction since the damping parameter a <C 1 . 
Although, the nanomagnet has small dimension along its thickness (i.e., I <^b < a and 
the demagnetization factors Ndx 3> Ndy > Ndz), magnetization cannot remain on the 
magnet’s plane {y-z plane, 0 = ±90°) since the interface coupling energy is a few orders 
of magnitude higher than the shape anisotropy energy. Thus, the magnetization keeps 
rotating in the 0 -direction, but it also traverses towards the anti-parallel direction in 
^-space (e.g., 9 ~ 180° to 0 ~ 0°) due to damping [see Eq. ([7])]. 

Note that exactly at 9 = 180° or 0°, the torque acting on the magnetization due to 
interface anisotropy [Eq. ([2])] is exactly zero, however, thermal fluctuations can scuttle 
the magnetization from these points to initiate switching. At the very start of switching, 
the initial orientation of magnetization is not a fixed value rather a distribution due to 
thermal agitations. Such distribution is considered during simulations. Hence, thermal 
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fluctuations affect the switching when magnetization starts to switch as well as during 
the course of switching. 

The energy dissipated in the nanomagnet due to Gilbert damping can be expressed 
as Ed = Pd(t)dt, where r is the switching delay and Pd(t) is the power dissipated at 
time t per unit volume given by 

Thermal held with mean zero does not cause any net energy dissipation but it causes 
variability in the energy dissipation by scuttling the trajectory of magnetization. 

3. Results 

As depicted in the Fig. [H the nanomagnets (Mi and M 2 layers) are made of Fe, while 
the ferroelectric P-layer is made of PbTiOz [IHl EH EEl EZ]. The spacer layer is made 
of Au and the thicknesses of the trilayer Mi/spacer/M 2 are 1/4/1 monolayers |18j . 
The Fe layer has a unit cell length of 0.287 nm and it has the following material 
parameters: saturation magnetization {Mg) = le6 A/m, and damping parameter (a) = 
0.01 [38l [39l 00]. The elliptical lateral cross-section has a dimension of 15 nm x 7 nm. 
The P-layer has a unit cell length of 0.388 nm and it has 5 layers in the vertical direction 
(x-direction, Fig. [1]) |T8] while the lateral dimensions are same as of the nanomagnets. 
The energy difference between the P-alignment and AP-alignment is 10 mev/atom [T8] . 
and the absolute value of energy turns out to be about 10 ev or 385 kT at room- 
temperature. This huge interface coupling makes potential landscape of Mi monostable 
at 6 = 180° or at 0 = 0° depending on the P-alignment or the AP-alignment in 
the trilayer, respectively. Hence, no spontaneous switching of magnetization between 
6 = 180° and 6 = 0° can take place. The shape anisotropic energy barrier is a few orders 
of magnitude lower than this interface coupling energy and thus consideration of shape 
anisotropy does not make any significant difference, however, it is included during the 
simulations. 

Modeling the P-layer as a parallel plate capacitor and using a relative dielectric 
constant of 1000 BH, the capacitance C of the layer becomes ~0.4 fF. If the P-layer is 
accessed with a 10 pm long silver wire with resistivity ~2.6 pH-cm [12], the resistance 
R becomes ~3 kfl. Therefore, the RC time constant is of the order of 1 ps. The 
ferroelectric PbTiO^ has a coercive voltage of 20 MV/m [13] and hence a voltage V ~ 40 
mv is required to switch the polarization. (Note that the voltage required to switch the 
traditional charge based devices is of the order of 1 V H.) Polarization switching is 
possible in less than 100 ps |15] and a voltage ramp with period T = 100 ps or more 
is considered to enforce the quasistatic (adiabatic) assumption (T ::§> RC). Without 
any adiabatic assumption, the metric CV^ is 0.5 aJ and hence the energy dissipation in 
the external circuitry that applies the voltage is miniscule. With 100 ps ramp period, 
the dissipation becomes a negligible value of 0.01 aJ [IHl [9]. Note that we do 

not calculate any standby leakage current through the thin ferroelectric since the device 
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Dynamics of 0 



Figure 2. A sample dynamics of magnetization while switching from 6 ~ 180° to 
0 ~ 0° in the presence of room-temperature (300 K) thermal fluctuations. The ramp 
period is 100 ps and the switching delay is 168.5 ps. The energy dissipation due to 
Gilbert damping is 1.42 aJ. 


operation is non-volatile, i.e., it is possible to tnrn off the voltage withont loosing the 
information. Also, dnring the active mode of operation, the tnnneling cnrrent is small 
(< 1 nA H), leading to negligible energy dissipation. Note that ferroelectric fatigne 
may make the coercive field higher [Ul HE], however, since the energy dissipation dne 
to applied voltage is miniscnle, it does not appear to be a bottleneck provided the 
polarization switches and the interface conpling between the P-layer and the trilayer 
still persists. 

Figure [2] shows a sample dynamics of magnetization in the presence of room- 
temperature (300 K) thermal fluctuations. The ramp period is considered to be 100 
ps and the switching has completed in less than 175 ps. During the course of switching 
magnetization has temporarily backtracked due to random thermal kicks, however, the 
strong interface anisotropy has enforced magnetization to switch from 6 ~ 180° to 
9 ~ 0°. 

Figure [3] plots the average switching delay versus average energy dissipation for 
different ramp periods (0.1 ns - 1 ns). A moderately large number of simulations 
(10000) in the presence of room-temperature (300 K) thermal fluctuations are performed 
to generate each point in the curve. When the magnetization reaches 6 < 5°, the 
switching is deemed to have completed. Note that as we vary the ramp period of 
applied voltage across the heterostructure slower, the switching also gets slower and less 
energy is dissipated in the switching process, elucidating the delay-energy trade-off for 
the device. The results show that switching in sub-nanosecond delay is plausible while 
dissipating energy of only ~1 aJ. The energy dissipation is a couple of orders of 

magnitude lower than the energy dissipation due to Gilbert damping and it decreases 
with the increase of ramp period since the switching becomes more adiabatic. The 
standard deviation in switching delay for ramp period of 0.1 ns is about 22 ps and it 
increases about twice when the ramp period is increased to 1 ns. At higher ramp period. 
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Figure 3. Switching delay-energy trade-off as a function of ramp period (upper axis). 
For a faster ramp, the switching becomes faster but the energy dissipation goes higher. 
Each point is generated from 10000 simulations in the presence of room-temperature 
(300 K) thermal fluctuations and the average values of switching delays and energy 
dissipations are plotted. For 0.1 ns ramp period, the average (max) switching delay is 
175.3 ps (330.8 ps), while the mean energy dissipation is 1.56 aJ. For a slower ramp 
with period 1 ns, the average (max) switching delay is 775.2 ps (1003.5 ps), while the 
mean energy dissipation is 0.58 aJ. 


thermal fluctuations have more time to scuttle the magnetization and cause variability 
in switching time. Simulations have been performed at an elevated temperature (400 
K) and the performance metrics switching delay and energy dissipation turn out to be 
similar (within 5%) compared to that of room-temperature (300 K) case [4^ . 

4. Discussions 

The giant magnetoresistance (GMR) [50l [51] of the trilayer in this interface-coupled 
structure is calculated to be of the order of 30% |T8l [52], which provides an way to 
read the magnetization states (P-alignment or AP-alignment). Although this GMR 
is not that high compared to tunneling magnetoresistance (TMR) [SS] [Ml ESI Ml EZ] , 
suitable design strategies can be possibly be devised to work with this moderate value of 
GMR and also it may be possible to increase the GMR by suitable material choice and 
design. It is also argued that even with the variance in the smaller thicknesses of the 
layers, it is still possible to couple the polarization and magnetization interracially in 
the proposed structures |T8] . The modeling of interface anisotropy is not limited to the 
way we performed here, however, any strong interface-coupled system would facilitate 
switching of magnetization from one state to the another. 

Due to the small lateral dimensions of these interface-coupled multiferroic 
heterostructures, it is possible to cram an enormous amount of devices on a single 
chip. Using an area density of 10“^^ cm“^, the dissipated power would be 10 mW/cm^ 
considering 1 aJ energy dissipation in a single nanomagnet with 1 ns switching delay 
and 10% switching activity (i.e., 10% of the magnets switch at a given time). Such 
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extremely dense and ultra-low-energy non-volatile computing systems can be powered 
by energy harvesting systems [58l [59l IMl 1 ^ - 

5. Conclusions 

In conclusion, we have performed an analysis over the switching dynamics of 
magnetization in interface-coupled multiferroic magnetoelectrics. We have modeled the 
interface anisotropy and calculated the performance metrics e.g., switching delay and 
energy dissipation as a function ramp period to elucidate the delay-energy trade-off. 
The results show that switching can take place in sub-nanosecond switching delay while 
expending ~1 aJ of energy. Also, strong interface coupling facilitates error-resiliency 
during the switching and allows to scale down the lateral area to very small dimensions. 
Due to these superior performance characteristics of interface-coupled multiferroics, it 
would be of immense interest to work out different possible theoretical designs and 
experimental implementations. Successful experimental implementations must tackle 
the issue of process variation at low dimensions. Processors built on such platform 
can harbinger unprecedented applications that can work by harvesting energy from the 
environment e.g., medically implanted device to warn an impending epileptic seizure by 
monitoring the brain signals, wearable computers powered by body movements etc. 
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SI. Switching delay-energy at T = 400 K 

To get an understanding on the performance metrics switching delay and energy 
at elevated temperatures, we have performed simulations at T = 400 K, which is 
particularly important for burn-in testing of the devices. Fig. [ST] shows the plot for 
switching delay versus energy dissipation at T = 400 K for different ramp periods (0.1 
ns - 1 ns). A moderately large number of simulations (10000) in the presence of thermal 
fluctuations (T = 400 K) have been performed to generate each point in the curve. When 
the magnetization reaches 6 < 5°, the switching is deemed to have completed. With the 
increase in switching delay, less energy is dissipated in the switching process, elucidating 
the delay-energy trade-off for the device. 


Ramp period (ns) 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure SI. Switching delay-energy trade-off as a function of ramp period (upper 
axis) at T = 400 K. For a faster ramp, the switching becomes faster but the energy 
dissipation goes higher. Each point is generated from 10000 simulations in the presence 
of thermal fluctuations (400 K) and the average values of switching delays and energy 
dissipations are plotted. For 0.1 ns ramp period, the mean (max) switching delay is 
172.3 ps (325 ps), while the mean energy dissipation is 1.56 aJ. For a slower ramp 
with period 1 ns, the mean (max) switching delay is 768 ps (970.2 ps), while the mean 
energy dissipation is 0.58 aJ. 


S2. Comparison of performance metrics for T = 300 K and T = 400 K 

We perform a comparison of switching delay (both mean and standard deviation) for 
the cases T = 300 K and T = 400 K. Figure [S2] depicts that the mean switching delay 
at a higher temperature T = 400 K decreases compared to the case at T = 300 K. This 
can be conceived by the reasoning that the initial deflection of magnetization due to 
thermal fluctuations increases at a higher temperature. Hence, magnetization is likely 
to start more far from the easy axis at a higher temperature for different trajectories, 
leading to the decrease in switching delay. However, this decrease in mean switching 
delay at T = 400 K is very small (less than 2%) compared to that of T = 300 K. 

Figure [S3] shows the standard deviation in switching delay for temperatures 300 
K and 400 K. The trend of standard deviation in switching delay with the increase in 
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Ramp vs mean switching delay 



Figure S2. Mean switching delay versus ramp period for T=300 K and T=400 K. 
Each point is generated from 10000 simulations in the presence of thermal fluctuations. 
At the higher temperature of T=400 K, the mean switching delay decreases compared 
to that of T=300 K. For 0.1 ns ramp period, the mean switching delay at T = 300 K 
(400 K) is 175.3 ps (172.3 ps). For a slower ramp with period 1 ns, the mean switching 
delay at T = 300 K (400 K) is 775.2 ps (768 ps). 


Ramp vs standard deviation in delay 



Figure S3. Standard deviation in switching delay versus ramp period for T=300 
K and T=400 K. Each point is generated from 10000 simulations in the presence of 
thermal fluctuations. At the higher temperature of T=400 K, the standard deviation 
in switching delay increases compared to that of T=300 K. For 0.1 ns ramp period, 
the standard deviation in switching delay at T = 300 K (400 K) is 21.8 ps (22.7 ps). 
For a slower ramp with period 1 ns, the standard deviation in switching delay at T = 
300 K (400 K) is 43.9 ps (45.4 ps). 


temperature shows an opposite trend to that of the mean. The standard deviation at 
higher temperature increases with temperature, which bodes well with the Eq. (3) in 
the main paper and also physically conceivable. This increase in standard deviation in 
switching delay at T = 400 K is quite small (less than 5%) compared to that of T = 
300 K. 

The mean energy dissipation decreases with increasing temperature and this 
decrease at T = 400 K is quite small (less than 3%) compared to the case of T = 
300 K. 















